Hamiltonian formalism in Friedmann cosmology and its quantization 
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We propose a Hamiltonian formalism for a generalized Friedmann-Roberson- Walker cosmology 
model in the presence of both a variable equation of state (EOS) parameter w(a) and a variable 
cosmological constant A(a), where a is the scale factor. This Hamiltonian system containing 1 
degree of freedom and without constraint, gives Friedmann equations as the equation of motion, 
which describes a mechanical system with a variable mass object moving in a potential field. After 
an appropriate transformation of the scale factor, this system can be further simplified to an object 
with constant mass moving in an effective potential field. In this framework, the A cold dark matter 
model as the current standard model of cosmology corresponds to a harmonic oscillator. We further 
generalize this formalism to take into account the bulk viscosity and other cases. The Hamiltonian 
can be quantized straightforwardly, but this is different from the approach of the Wheeler-DeWitt 
equation in quantum cosmology. 

PACS numbers: 98.80. Jk,45.20.Jj, 03. 50.-Z 
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I. INTRODUCTION 



Since the current accelerating expansion of our Uni- 
verse was discovered [l[ around 1998 and 1999, theoret- 
ical physicists have devoted increasingly more attention 
to the Friedmann-Roberson- Walker (FRW) model as a 
standard framework in cosmology study. The A cold dark 
matter (ACDM) model as the standard model of cosmol- 
ogy so far fits well with observational data whereas it 
has had some serious theoretical problems. To make a 
comparison to the ACDM model, physicists have built 
many cosmological models that are able to give out the 
effective Friedmann equations with variable cosmological 
constant. To quantize the Friedmann equations, the com- 
monly used theory is the Wheeler-DeWitt equation Q, 
which has been studied and applied widely in quantum 
cosmology Q. Starting from the Hilbert-Einstein action 
with the Roberson- Walker (RW) metric, the correspond- 
ing Hamiltonian TL can be obtained. Then the Friedmann 
equation plays the role as the constraint Ti = 0, which 
leads to the Wheeler-DeWitt equation. In the present 
work, we consider the Friedmann equations as basic equa- 
tions and find a Hamiltonian system that gives Fried- 
mann equations as classical equations of motion without 
constraint. 

Many Ansatze of the variable cosmological constant 
have been studied in the literature, such as Refs. 0, d, 
0, 0|- Moreover, some models motivated from the string 
theory give an effective cosmological term when reduced 
to the FRW framework. We assume that the equation 
of state (EOS) parameter w = p/p can also be variable, 
which means that the contents of the Universe, except 
the cosmological term, are generalized to a nonperfect 
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fluid, or perfect fluid as a special case. In observational 
cosmology, the redshift z is regarded as an observable 
quantity and related to the scale factor a by z = ao/a—1. 
Therefore, we investigate a general case that both the 
EOS parameter w and the cosmological constant A can 
be functions of the scale factor a, and take into account 
the bulk viscosity. 

As an extension of the problem, we construct a Hamil- 
tonian formalism for a system described by the following 
equation: 

q = h{q)q 2 +vq + h(q), 

where fi{q) and /2(<?) are arbitrary functions, and r\ is 
constant. Also it can be regarded as a generalization 
of the damping harmonic oscillator. The corresponding 
Hamiltonian describes an object with variable mass mov- 
ing in a potential field. After an appropriate canonical 
transformation, this system can be further simplified to 
an object with constant mass moving in an effective po- 
tential held. Thus, differential models in the FRW frame- 
work are characterized by their effective potentials. This 
is a general formalism and it can be applied to many cos- 
mological models, for example, that the ACDM model 
corresponds to a harmonic oscillator. Since the quanti- 
zation of Friedmann equations can provide an insight to 
quantum cosmology as a glimpse of quantum gravity, we 
also make some remarks on the quantum case, which pro- 
vides a correspondence between cosmology and quantum 
mechanics. 

The paper is organized as follows. In Sec. II we present 
a generalized FRW model and the corresponding Hamil- 
tonian to describe the Friedmann equations. Then we 
find a canonical transformation to further simplify the 
problem, and give some examples and special cases. In 
Sec. Ill we show that our framework can also be applied 
in the dissipative case with bulk viscosity. In Sec. IV 
we turn our attention to the relation to the observable 
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quantities and review some issues of the Bianchi identity. 
In Sec. V we make some remarks on quantum cosmology 
from our approach. In the last section we present the 
conclusion and discuss some future subjects. 



II. HAMILTONIAN FORMALISM 

A. Hamiltonian description of the Friedmann 
equations 

We consider the RW metric in the flat space geometry 
(fc=0) as the case favored by current cosmic observational 
data: 



ds 2 



-dt 2 + a{t) 2 {dr 2 +r 2 dQ 2 ), 



(1) 



where a(t) is the scale factor. The energy-momentum 
tensor for the cosmic fluid can be written as 



Tpv = {p + P)U^U V + {p + pa) g,i V , 



(2) 



where pa = A/(8irG) is the energy density of the cos- 
mological constant. Thus, Einstein's equation R^ v — 
\g^R — SttGT^ contains two independent equations: 



8ttG 



A 

3 : 



3 

4nG. n . A 
-— - (p + 3p) + -. 



(3a) 
(3b) 



The EOS of the matter (cosmic fluid except the cosmo- 
logical constant) is commonly assumed to be 



P = (7- !)/>• 



(4) 



Cosmologists usually call Eq. (|3"a|) as the Friedmann 
equation and Eq. (|3b[) as the acceleration equation in the 
literature, whereas for simplicity we name both Eqs. (|3ap 
and (|3b[) Friedmann equations here. For generality, we 
assume that both 7 and A are functions of the scale factor 
a, thus we call it the generalized FRW model. Combining 
the Friedmann equations with the EOS, we obtain 



37(a) -2 a 2 7(a)A(o) 



(5) 



which determines the evolution of the scale factor. 

We regard Eq. §5§ as a basic starting point; therefore, if 
the dynamical equation for the scale factor can be written 
as that form, the present framework can be valid. If 
the Newton constant G is constant and the cosmological 
constant A is variable, the energy-momentum tensor for 
the matter cannot individually conserved 0, 0], which 
implies an interaction between the matter and vacuum 
energy. In the following, we assume G to be constant 
until Sec. IV. 



Our aim is to find a Hamiltonian description of Eq. ([5]) 
as the classical equation of motion. We start from the 
following Lagrangian 



C(q,q) = ~M(q)q 2 -V(q), 
and the corresponding Hamiltonian thus is 



n{ q , P ) 



p 



2M(q) 



V{q), 



(6) 



(7) 



with the canonical Poisson bracket {q,p} = 1. One can 
check that the equation of motion for Eq. |6j) or |7|) is 



q _ 1 d In M q 2 
q ~ ~2 dlnq f ~ ~Wq 



1 dV 



(8) 



This equation possesses the same form as Eq. ([5]). There- 
fore, by comparing Eq. ([5]) with Eq. (jHJ), we can take a as 
the general coordinate and solve the functions M(a) and 
V(a). Then the Lagrangian C = \M{a)q 2 — V(a) with 

M = exp (^J 3l ~ 2 da \ , v = -\ J M~fAada, (9) 

gives Eq. ([5|) as the equation of motion. For some spec- 
ified functions 7 = 7(a) and A = A (a), the above in- 
tegrations can be evaluated out to give M(a) and V(a) 
explicitly. 

Now we can see that the generalized FRW model essen- 
tially corresponds to an object with variable mass M(a) 
moving in a potential field V(a). In the following, we 
will show that this picture can be further simplified as 
an object with constant mass moving in an effective po- 
tential field V{4>), after an appropriate transformation of 
the scale factor. 



B. Canonical transformation 

The above problem can be generalized as the Hamilto- 
nian description of the nonlinear equation 



fi(q)q 2 + /2(g), 



(10) 



where fi(q) and /2(a) are two specified functions. 
This equation can be derived by the Lagrangian C = 
\M{q)q 2 - V(q) with 

M = exp (-2 J h(q)dq) , V = -J Mf 2 (q)dq. (11) 

We define a new variable 4> as (see Appendix) 

cf> = J exp(-J h{q)dq^dq. (12) 

This transformation can eliminate the q 2 term and gives 
the equation for the variable <p as in 



4> = fziq) exp - / fx(q)dq 



(13) 
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where q — > <f> denotes using Eq. JT2J) to change the vari- 
able q to 4>. Since there is no (j> 2 term in Eq. (|13[) . this can 
be regarded as a partial linearization. Therefore, the sys- 
tem of Eq. (fT0|) transformed to Eq. (fl~3|) can be described 
by the Lagrangian 



Here the parameters 7, A, and m are all constants. This 
equation possesses the same form of Eq. (|10[) . By defining 
4> as Eq. (fT2f and using Eq. (fT5| . we obtain the effective 
potential as 



with the potential as 



f 2 {q) exp 
/2(g) exp I -2 



fi(q)dq 



(14) 
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7 2 A^ 2 



A: 



37 



_ /3t0\ 
m \ 2 / 



2-2m/37 



(21) 



fi(g)dq I dg 



.(15) 



The simplification of the problem by Eq. (fl"2"|) is essen- 
tially the canonical transformation 



for 7^0 and m 7^ 37. 

Another example is the Friedmann equations during 
the inflation era. In the study of inflation, we usually 
use the conformal time r instead of the comoving cosmic 
time t. Here we assume that a constant term —po is in 
the EOS during inflation. Thus the Friedmann equations 
combined with the EOS p = —p — po yield 



a" a 12 k 2 

7? = 2 ^ + T Po ' 



(22) 



where the prime denotes a derivative with respect to r, 
q^(/>, Pq^Pc/,, H(q,Pq) -> U{(j>,p <t> ), (16) and n 2 = 8ttG. By defining cj) = -1/a, the equation for 

<j) is 4>"4> + (n 2 /2)po = 0. The effective potential is thus 



where p q = M(q)q, p^ = <fi, and H(^,p^,) = \p\ + V(<j>). 
Therefore, the classical and quantum properties of differ- 
ent models are characterized by the effective potentials. 

For Eq. ([5]) as a special case, the new variable 4> is 
given by 



exp 



37- 2 
2a 



da da. 



(17) 



V(<P) 



-p ln\(j>\. 



(23) 



Moreover, if we add the curvature term in this case, 
corresponds to a cj> 2 potential. 

III. BULK VISCOSITY 



C. Some examples 

We will give some special cases of the above general 
framework to show some applications. If both 7 and A 
are constant for a simple case, the integrations in Eq. (117|) 
can be evaluated out as 



37 

In a, 7 = 0. 



(18a) 
(18b) 



Now we consider 7^0 for example. The special case 7 = 
1 corresponds to the ACDM model. The equation for <fi 
can be obtained as <j>— §7 2 A(/> = 0, and the corresponding 
Lagrangian is 



C 



^ 2 + l^ 2 . 



(19) 



We can see that the simplest model in cosmology just 
corresponds to a harmonic oscillator after linearization. 
In particular, this is a upside-down harmonic oscillator 
for the asymptotic de Sitter Universe. 

We can add the curvature effect to the ACDM model, 
which is described by the special case m = 2 of the fol- 
lowing equation: 



ii _ 37 - 2 a 2 7A 
a ~ 2 a? + ~2~ 



(20) 



We assume that the cosmic fluid possesses some dis- 
sipation effects. Since the sheer tensor = for RW 
metric, the sheer viscosity does not contribute to the evo- 
lution in Friedmann cosmology. The energy-momentum 
tensor for nonperfect fluid concerning bulk viscosity in 
the right-hand side of Einstein's equation is given by d,Q 



Tpu = pU^Uy + (p - Co9)h^, 



(24) 



where h^ v = + U^U U is the projection operator, 
9 = J7." = 3d I a is the scalar expansion, and £ is the 
bulk viscosity coefficient. Consequently, Eq. ([5]) should 
be modified as 



a 
a 



- 37(a) " 2 ^ + 12 7 rGCo^ + 7(a)A(a) . (25) 

n ^ n 9 



where both 7 and A can be functions of a for generality, 
and Co is constant. We also find a Hamiltonian 



H(a,p a ,t) 



2M(a,t) 



+ V(a,t), 



(26) 



with the Poisson bracket {a,p a } = 1 to give Eq. (j2"5|) as 
the classical equation of motion. The functions in this 
Hamiltonian are given by 



M 



exp 



3 7 -2 



da - l2irG( Q t 



V = -- I M-fAada. 



(27a) 
(27b) 
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Although a dissipative system cannot be described by 
a conservative Hamiltonian generally, one can directly 
check that the classical equation of motion for the Hamil- 
tonian Eq. (f2"6"| is Eq. (j2"5)) . As a special case, the equa- 
tion for a damping harmonic oscillator can be derived by 
the Caldirora-Kani (CK) Hamiltonian fiot ] - 

The above problem can be generalized to construct a 
Hamiltonian system for the equation 



fi(q)q 2 +W + / 2 (9), 



(28) 



where r\ is constant. It can be derived by the Hamiltonian 
H(q,p,t) = iM(g,t)-V + %«) with 

M = exp (-2 f h{q)dq - rjt] , V = - [ Mf 2 (q)dq. 



Similarly, by using the new variable <f> defined by 
Eq. (fT2|) , the equation for <f> is 



(j) = V(j)+ f 2 (q) exp - / fi(q)dq 



(30) 



9^0 



Now we consider a very special case that both 7 and A 
are constant; then cj> defined by Eq. (|18ap satisfies 



4> - 12ttGCo^ - ^7 2 A0 = 0, 

which describes a damping harmonic oscillator. 
The damping harmonic oscillator 

8V(q) 



Mq = -m(t) 



()q 



(31) 



(32) 



has been studied in quantum mechanics. The CK Hamil- 
tonian 



H 



1 

2M' 



-r)t/M p 2 



^M^e^q 2 , 



(33) 



with the commutation relation [q,p\ = ih, can yield the 
dissipation equation (|32p through the Heisenberg equa- 
tion [10| | . Our work can be regarded as a generalization 
to the case of variable mass. It is the variable mass that 
generates a nonlinear term in the equation of motion that 
describes the generalized FRW model. 

In our previous work [9j], we have proposed an EOS as 

2 _ 2 



P= (7- 



(34) 



where the parameters 7, Ti and T 2 are constants. Com- 
bining the Fricdmann equations with this more practical 
EOS, we obtain the dynamical evolution equation for the 
scale factor as 

a 37 - 2 a 2 la 1 

- = ~o 2+^F- + Tf2- 35 

a 2 or T\ a 1 2 

This model possesses a large variety of properties, such as 
that we have found a scalar field model which is equiva- 
lent to the above EOS. For related works on the modified 
EOS, see Ref. @, OH, El d • The present work can also 
be regarded as a generalization of the EOS to 7 = 7(a) 
and T 2 = T 2 (a). And the corresponding Hamiltonian 
formalism for this system can be constructed similarly. 



IV. RELATIONS TO THE OBSERVABLE 
QUANTITIES 

The observations of the supernovae (SNe) la have 
provided the direct evidence for the cosmic accelerat- 
ing expansion of our current Universe 0. A bridge 
between the cosmological theory and the observation 
data is the H-z relation, where H = a/a is the Hub- 
ble parameter and z is the redshift. For example, the 
ACDM model in cosmology can be described mainly as 
H 2 (z) — H 2 [Cl m (l + z) 3 + 1 - Cl m ], where Q m is the mat- 
ter energy density. This model fits the observational data 
well and provides the cosmological constant as the sim- 
plest candidate for dark energy. In a sense, the different 
cosmological models are characterized by the correspond- 
ing H-z relations. 

There is also a systematic way to construct the Hamil- 
tonian starting from the general model 



H 2 



/(«), 



(36) 



where /(a) is a specific function of the scale factor a, 
according to the model. By differentiating Eq. (|3"6"|) , we 
obtain that it is a solution of the following equation: 



3 7 -2 a 2 , 37/(0) , a/'(a) 



(37) 



which possesses the same form of Eq. or (jTU]) . The 
corresponding coefficients are given by 



Mo) 



37 



2 a 



Ma) 



3 7 a/ (a) , a 2 f'(a) 



(38) 



Then by applying Eq. (|lip we can obtain the correspond- 
ing Hamiltonian. Therefore, even if the EOS for a cosmo- 
logical model is not explicitly linear in p, the Hamiltonian 
formalism in the present work can also be applied if the 
effective Fricdmann equation H 2 = f(a) can be given out 
for that model. 

Many approaches such as modified gravity [l4| can 
be reduced to effective Friedmann equations in the form 
H 2 = f(a). Since ACDM model fits the SNe la data well, 
the reasonable cosmological models should be reduced to 
Friedmann cosmology in an effective way and give out the 
right H-z relation, in order to make a comparison with 
the ACDM model. In our case, the Friedmann equations 
in terms of the Hubble parameter can be written as 



aH- 



,dH 
da 



3 -1h 2 



A(a) 



(39) 



Here 7 is assumed to be constant for simplicity. This 
equation is linear in H 2 and the effective term A(a) is 
an inhomogeneous term. The solution in terms of H (z) 
concerning the initial condition H (0) = Hq is given out 
by 



H(z) 2 ^H 2 (l + zf 



1 



A(z')(l + z)-^- x dz' 



(40) 



5 



In the power-law ACDM model, the contributions of dif- 
ferent components are separated in H 2 , such as a con- 
stant for the cosmological constant, and a (1 + z) 2 fac- 
tor for the curvature term. But in the general case, the 
contribution of the matter cannot be separated from the 
above solution. This problem is related to the conserva- 
tion law of the matter, which has been investigated in 
Refs. [Si]. 

The Bianchi identity for the energy-momentum tensor 
Eq. © gives 



PA + p + 3H(p + p) = 0, 



(41) 



which implies that energy transfer will exist between the 
matter and the vacuum energy. An intuitive idea has 
been proposed that if both G and A are variable, the 
ordinary energy-momentum tensor can be individually 
conserved, i.e., p + 3H(p + p) = @. This is achieved by 
combining the Bianchi identity for the variable G and A 
model 



-[G(p A + p)] + 3GH(p + p) = 0, (42) 



with the following constraint: 



{p + PA )G + Gp A = 0. 



(43) 



The authors of Ref. |6j assume that both the Newton con- 
stant G and the cosmological constant A are functions of 
a scale parameter p and apply the renormalization group 
approach to cosmology. If G{p) evolves by a logarithmic 
law and pa(p) evolves quadratically with p, then this 
picture can explain the evolution of the Universe, and at 
the same time, the variable G can explain the flat rota- 
tion curves of the galaxies without introducing the dark 
matter hypothesis. 



V. REMARKS ON QUANTUM COSMOLOGY 

We have obtained a classical Hamiltonian formalism 
of the Friedmann equations. Generally, once a Hamilto- 
nian is obtained, the system can be quantized straightfor- 
wardly by replacing the Poisson bracket with the commu- 
tation relation [q,p] = i. However, we need to take into 
account the ambiguity in the ordering of noncommuting 
operators q and p. For simplicity, we ignore the order- 
ing ambiguity here. In terms of the new variable 0, the 
corresponding Schrodinger's equation can be written as 



(44) 



where p^, = —id$. To make a comparison between our ap- 
proach and the Wheeler-DeWitt equation, we only take 
the ACDM model as a very special case for an illustrative 
example. The corresponding Hamiltonian for Eq. Q19p in 
the case 7 = 1 is 



2a 



1 



:Aa 3 , 



(45) 



where p — ah. In the approach of the Wheeler-DeWitt 
equation, TC = is a constraint [2l Ha|. thus the quanti- 
zation gives (d 2 + ^a 4 ) v I / (a) = 0. This is an anharmonic 
oscillator with zero energy eigenvalue. In our case, the 
Hamiltonian is nonzero and proportional to the matter 
energy density, which we show in the following. The so- 
lution of Eq. (09j) with A(a) = A/2 is 



H 2 



A 



Hq [£! TO a 3 + 1 — f2 m ], 



(46) 



where ft m = 1 — A/(3_ffg). Therefore, the Hamiltonian 
can be calculated as 



n 



a -( a - 2 -±) = l -H 2 n m . 

2 V a 2 3 J 2 



(47) 



After a canonical transformation by Eq. ([T5)| . the 
Schrodinger's equation in terms of <p becomes 



YaW 



-A<f> 2 



*(0) = E^((j)). 



(48) 



Thus, for the asymptotic de Sitter Universe, the ACDM 
model corresponds to an upside-down harmonic oscillator 
in our formalism. Such an oscillator also appears in the 
matrix description of de Sitter gravity (lrij . 

We can transform the de Sitter Universe to the dual 
anti-de Sitter Universe by employing the scale factor du- 
ality [l7j . which has been found that a — > a -1 gives 
H — > — H and other consequences. The duality for 
Eq. ([5]) is given by 



a — > a 



7 



-7, A 



-A, 



(49) 



It can be checked easily that Eq. ([5]) is invariant under 
these transformations. If we use the dual scale factor a -1 , 
the corresponding potential becomes V(<f>) — +|7 2 A^ 2 . 
In fact, quantization in de Sitter spacetime is one of the 
major difficulties of string theory at one time (though 
this picture has changed a little bit after Kachru-Kallosh- 
Linde- Trivedi theory appeared). It seems that quantiz- 
ing de Sitter cosmology is no difference, since the time 
variable used is the same, and it is known that there is no 
global timelike coordinates in de Sitter spacetime. Some 
quantum effects of a scalar field in de Sitter background 
can be found in Ref. [181 ]. 



VI. CONCLUSION AND DISCUSSION 

We have proposed a systematic scheme to describe the 
Friedmann equations through a Hamiltonian formalism. 
The generalized FRW model accompanied by both vari- 
able EOS parameter and variable cosmological constant 
admits a Hamiltonian description without constraint. Af- 
ter an appropriate canonical transformation, the system 
can be significantly simplified to an object moving in 
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an effective potential field. The bulk viscosity can also 
be taken into account by a time-dependent Hamiltonian. 
Some examples are given explicitly such as the ACDM 
model, the curvature term effect, and the inflation period. 
The quantization of the system provides a new approach 
to study the potential quantum cosmology which is an 
intriguing topic in theoretical physics research. 

We shall discuss some possible future developments of 
our work. As we have claimed, the formalism in this work 
can be applied to a large variety of cosmological models. 
By solving the Schrodinger equation Ti.((j),p ! f,)'i l = Ety, 
the cosmological wave function can be obtained for a 
specific model. Here we consider the curvature effect, 
for example, which is described by the potential Eq. (j2"Tj) 
with parameters A = 0, 7 = 1, and m = 2. The cor- 
responding Schrodinger equation can be solved in terms 
of the biconfluent Heun equation (BHE) 19]. We can 
also start from the effective Lagrangian and study the 
observational effects when we modify the potential. We 
believe that our formalism would give a new perspective 
to the potential study of quantum cosmology physics. 
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APPENDIX A: MATHEMATICAL NOTES 

A more general correspondence between a Hamiltonian 
and its equation of motion is given in Ref. 19]. The 
equation of motion of the Hamiltonian 



is given out by 



Q = 



ld\nP q 2 (d\nf <91nP 



2 dq q 2 V 9t dt 
\dq2P T dtP 8q 



(A2) 



In the mathematical aspect, Eq. (|28|) can be further gen- 
eralized to the following equation: 

q = F 1 (q, t)q 2 + F 2 (q, t)q + F 3 (q, t), (A3) 

however, here the coefficients F\ and Fi are not com- 
pletely independent. Comparing with Eq. (|A2[) , we can 
see that the condition 2d t Fi(q,t) = 9 g F 2 (g,t) must be 
satisfied for consistency. In the present work, both f\ (q) 
and 77 have safely satisfied this condition. 

We shall explain why we choose the transformation as 
in Eq. (fT2"|) . Starting from the following equation 



q = h(q)q 2 + m + /2(g), 



(A4) 



we expect that after an appropriate change of variable 
<p(q), the above equation can be transformed as 



4> = r}()> + g{4>). 



(A5) 



By differentiating </>((/), we obtain cj) = (j)'q, and <fi = 
4>" q + <p'q\ where the prime denotes a derivative with 



respect to q. Substituting <fi, (f>, and <fi into Eq. (|A5[) , we 
obtain 



> .2 , ■ , 9 

J7« 



(A6) 



1 



H(q,p,t) = j^- (Po (q, t)p 2 + P 1 (q,t)p + P 2 (q,t)) 



(Al) 



Now it turns out that by defining —(/>"/$' = f\(q), which 
can be solved as the form Eq. (fl2")) . the q 2 term can be 
eliminated. 
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